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Henry’s Law solubility constants for HC1 have been measured for liquid 
H 2 S0 4 -HN0 3 -HC1-H 2 0 solutions; the results are in good agreement with 
predictions from published semiempirical models. The C10N0 2 + HC1 
reaction on the surfaces of such solutions with compositions simulating 
those of stratospheric aerosols has been investigated; as the composition 
changes following the temperature drop characteristic of the high-latitude 
stratosphere the reaction probability y increases rapidly. Furthermore, the y 
values remain essentially unchanged when HN0 3 uptake is neglected; the 
controlling factor appears to be the solubility of HC1. These results corrobo- 
rate our earlier suggestion that supercooled liquid sulfate aerosols promote 
chlorine activation at low temperatures as efficiently as solid polar strato- 
spheric cloud particles. 


The importance of heterogeneous reactions involving C10N0 2 , HC1, N 2 O s and H 2 0 
on the surfaces of polar stratospheric cloud (PSC) particles has been well established in 
connection with the observation of large-scale ozone depletion during the Antarctic 
springtime. 1,2 Laboratory studies have shown that both solid and liquid particles can 
promote heterogeneous reactions such as the following, converting the relatively inert 
chlorine reservoir molecules (C10N0 2 and HC1) into species (Cl 2 , HOC1) which are 
readily photolysed to yield chlorine radicals: 

C10N0 2 + H 2 0 - HOC1 4- HN0 3 (1) 

HOC1 + HC1 Cl 2 + H 2 0 (2) 

C10N0 2 + HC1 - Cl 2 + HN0 3 (3) 

In particular, heterogeneous reactions on sulfuric acid aerosols have recently been under 
intense study since these particles may play a role in mid-latitude ozone depletion. 3-6 
Hanson et al 1 have described a framework for applying the results of these experimental 
studies to the modelling of the actual micrometre sized liquid sulfuric acid aerosols 
found in the stratosphere. It was concluded in these studies that reactions (l)-(3) become 
important at sulfate aerosol compositions corresponding to low atmospheric tem- 
peratures ( T < 200 K), in agreement with our earlier suggestion that stratospheric 

aerosol particles need not freeze in order for chlorine activation reactions to proceed 
efficiently at high latitudes. 8 

It is important to determine accurately the composition of the aerosols as a function 
of the temperature and partial pressures of their constituents in order to model the 
various possible heterogeneous reactions adequately. Several groups have investigated 
uptake of HN0 3 9 and HC1 4, 10,11 by sulfuric acid solutions. Vapour pressures of these 
constituents dissolved in bulk solutions have also been measured. 12 In addition, some of 
this data has been used to construct semi-empirical thermodynamic models which are 
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intended to predict sulfate aerosol composition. 1314 Although it has been shown that 
HN0 3 and HC1 uptake by stratospheric sulfate aerosols is negligible for equilibrium 
temperatures above 210 K, the solubility of these species increases rapidly with decreas- 
ing temperature. Field observations of total aerosol volume at high northern latitudes 
have commonly been interpreted as evidence for solid PSC formation, since the aerosol 
volume is observed to rise quickly over a temperature range near the nitric acid tri- 
hydrate (NAT) frost point. 15 However, these observations are more consistent with the 
existence of supercooled liquid sulfate aerosols with substantial nitric acid content than 
with the existence of solid PSC particles. 16 

Since the reaction probability for the C10N0 2 + HC1 reaction depends very sensiti- 
vely on the HC1 concentration in the aerosol, it is particularly important to determine 
accurately the solubility of HC1 in solutions with compositions of atmospheric relevance. 
The existing results for the HC1 solubility in sulfuric acid solutions are in conflict; in 
particular, the Henry’s Law coefficients determined earlier in our laboratory 12 deviate 
significantly from similarly extensive measurements of Hanson and Ravishankara 4 and 
from the thermodynamic model predictions of Luo et a/. 13 

Here we describe further measurements designed to resolve the discrepancy in the 
solubility results for the ternary H 2 S0 4 -H 2 0-HC1 system, and report the first measure- 
ments of Henry’s Law coefficients for HC1 for the quaternary (H 2 S0 4 -HN0 3 - 
H 2 0-HC1) system. We also report measurements of C10N0 2 + HC1 reaction probabil- 
ities on the surface of such ternary and quaternary supercooled solutions with composi- 
tions relevant to the stratosphere. 


Experimental 

The HC1 vapour pressures and the reaction probabilities for the C10N0 2 + HC1 reac- 
tion were measured using a vertical wetted-wall flow tube coupled to a differentially 
pumped molecular-beam-sampling quadrupole mass spectrometer. The vertical flow 
tube is similar to that employed by Utter et al} 1 for measurements of ozone uptake on 
aqueous surfaces. 

The top of the vertically aligned flow tube (id 2.4 cm) was fitted with a liquid delivery 
system that contained inlets for the helium carrier gas, a movable injector and the sul- 
furic acid solution. The delivery system consisted of an annular cup which, when filled 
with liquid, overflows to create a falling cylindrical film of acid on the inner wall of the 
flow tube. This wetted-wall technique has been widely used in the field of chemical 
engineering, and Danckwerts 18 has described the fluid dynamics of such flow. In order 
to maintain cold, stable temperatures, it was necessary to shut off the acid solution flow 
during actual measurement of HC1 vapour pressures. The temperature of the acid was 
assumed to be the temperature of the coolant (which was sampled at the inlet to the flow 
tube) when the HC1 signal was observed to stabilize (of the order of a few minutes). The 
flow-tube wall was rinsed with fresh solution for each measurement in order to minimize 
the possible effects of dehydration : there is insignificant dehydration of such acid solu- 
tions on such short timescales, as shown by Hanson and Ravishankara. 19 The carrier 
gas was sampled at the end of the flow tube and directed to the detection system, while 
the slowly moving acid film was collected in a bulb maintained at temperatures below 
that of the flow tube in order to prevent backstreaming of the vapours. 

Most measurements were made with flow tube pressures of ca. 1 Torr and carrier gas 
flow velocities of ca. 1000 cm s _1 . The saturation of the carrier gas with HC1 vapour was 
checked by directing the helium flow through the moveable injector; the observed HC1 
signal did not change with injector position, thus confirming the equilibrium condition. 
In order to determine absolute values for the HC1 vapour pressures, the mass spectro- 
meter was calibrated using the extrapolated HC1-H 2 0 data of Fritz and Fuget. 20 The 
detection limit for HC1 (monitored at mass 36) using these methods was ca. 10“ 7 Torr. 
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C10N0 2 partial pressures of ca. 1(T 6 Torr were used for the kinetics measurements 
and the possible reactive depletion of HC1 was investigated for each solution. Because 
the electron-impact ionization scheme results in the production of N0 2 + for both 
C10N0 2 and HN0 3 , it was necessary to detect the reaction product, Cl 2 , for the quat- 
ernary solutions (since the ambient HN0 3 vapour pressures were high enough to inter- 
fere with the C10N0 2 signals). Reaction probabilities were determined using standard 
cylindrical flow tube analysis techniques. 21 

Bulk solutions were prepared by diluting 96% H 2 S0 4 , 70% HN0 3 and 37% HC1 
with deionized water. The solutions were spot-checked for composition with density and 
standard acid-base titration measurements. 


Results and Discusssion 


HC1 Vapour Pressures 

Henry’s law solubility constants for HC1 were determined for two ternary solutions and 
for two quaternary solutions over a range of temperatures (ca. 205-235 K) of strato- 
spheric relevance. A 50% H 2 S0 4 ternary solution was chosen for study to allow direct 
comparison with previous measurements. The other solutions were chosen in order to 
test the thermodynamic model of Luo et al. 13 for compositions of atmospheric interest. 
Neglecting uptake of HN0 3 , it may be inferred that at a stratospheric temperature of 
ca. 190 K at 50 mbar the aerosols will contain 43% H 2 S0 4 by weight, assuming 5 ppm 
by volume H 2 0. We therefore measured solubility constants for this composition, for 
which the uptake of HC1 is almost surely large enough for the C10N0 2 + HC1 reaction 
to be quite important. In order to spot-check the Luo et al. model further, we chose two 
quaternary solutions that correspond to aerosol compositions at equilibrium tem- 
peratures of 198 and 196 K at 100 mbar, with 10 ppb by volume HN0 3 and 2 ppb by 
volume HC1; these solutions contain significant levels of nitric acid (3.5% and 12.5%, 
respectively). 

Table 1 contains the measured Henry’s Law solubility constants ( H *) for the four 
solutions considered. In Fig. 1, the various results for the 50% H 2 S0 4 solution are 
presented. In earlier work in our laboratory 12 the vapour pressure data was para- 
metrized as a function of temperature, and it is shown in Fig. 1 as a dashed line. The 
results from the model of Luo et al. 13 are shown as a solid line, and each data point 
from Hanson and Ravishankara 4 Williams and Golden, 11 and this work is shown indi- 
vidually. Our new results form a consensus with two of the previous measurements and 
with the model predictions, while the Zhang et al. 12 data lies well outside the range of 
the consensus findings. The reason for the discrepancy is not clear; it is most likely a 
consequence of HC1 calibration errors in our earlier work. In Fig. 2, the H* values for 
the 43% H 2 S0 4 solution are plotted together with the model predictions of Luo et al.; 13 
as can be seen in the figure, the agreement is very good. 

In Fig. 3, our measured HC1 Henry’s Law constants for the two quaternary solutions 
are presented together with the predictions from the thermodynamic model; once again, 
for these atmospherically relevant compositions, the agreement is excellent. These mea- 
surements are not intended to be an exhaustive test of the accuracy of the model; never- 
theless, as has been shown to be the case for the ternary solutions it appears that the 
model also predicts the HC1 solubility for quaternary systems very accurately. 


Reaction Probabilities for HCI 4* C10N0 2 

The compositions of the bulk solutions used for the kinetics investigations simulated 
those of typical stratospheric aerosols. These compositions were determined by using the 
Luo et al. 13 model to search for mixtures whose vapour pressures correspond (within 
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Table 1 Experimental Henry’s law solubility constants 


composition 

77 K 

H */ mol dm 3 atm’ 1 

50% H 2 S0 4 

225 

4.0 x 10 5 

9.2 x 10 3 mol dm -3 HC1 

222 

5.2 x 10 5 


219 

7.2 x 10 5 


216 

9.8 x 10 5 


214 

1.2 x I0 6 


211 

1.8 x 10 6 


208 

2.3 x 10 6 

43% H 2 S0 4 

226 

2.5 x 10 6 

3.0 x 10 2 mol dm -3 HC1 

222 

4.1 x 10 6 


218 

6.8 x 10 6 


214 

9.2 x 10 6 


208 

1.8 x 10 7 

48% H 2 S0 4 

233 

8.8 x 10 4 

3.5% HN0 3 

231 

1.1 x 10 5 

3.9 x 10’ 3 mol dm 3 HC1 

226 

1.7 x 10 5 


222 

2.7 x 10 5 


216 

5.4 x 10 5 

36.2% H 2 S0 4 

228 

5.5 x 10 5 

12.5% HN0 3 

223 

1.2 x 10 6 

6.2 x 10’ 3 mol dm 3 HC1 

218 

1.7 x 10 6 


213 

2.8 x 10 6 


208 

4.7 x 10 6 


Tf K 

230 220 210 200 



Fig. 1 Henry’s Law constants as a function of temperature for a 50% H 2 S0 4 solution. ( — 

data from Zhang et a/., 12 and ( ) from the semi-empirical model of Luo et a/. 13 A, our data, 

data from Williams and Golden 11 and #, from Hanson and Ravishankara. 4 
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7/ K 

230 220 210 200 



Fig. 2 Henry’s Law constants as function of temperature for a 43% H 2 S0 4 solution. The line 
represents estimates from the Luo et al. model 13 and A, represent our data. 



4.2 4 4 4.6 4 6 

10 3 K/T 


Fig. 3 Henry’s law constants as function of temperature for the two quaternary solutions studied 

in this work. The lines represent estimates from the Luo el al. model 13 and A, our data. ( ) and 

▲, 36.2% H 2 S0 4 -12.5% HNO., solution, ( ) and A, 48.0% H 2 S0 4 3.5% HN0 3 solution. 
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10%) to various atmospheric equilibrium temperatures T(eq), assuming 5 ppm by 
volume H 2 0 and 2 ppb by volume HC1 at 100 mbar ( ca . 16 km altitude), with no HN0 3 
in Table 2 and 10 ppb by volume HN0 3 in Table 3. As the equilibrium temperature is 
lowered from 199 to 195 K, the weight fraction of HNO s in the solution increases from 
2.2% to 28.3%, and the amount of HC1 increases by more than a factor of 30. For 
aerosols at equilibrium with an atmosphere of 50 mbar total pressure and the same 
mixing ratios, similar compositional changes are expected at temperatures shifted by ca. 
4 K below the 100 mbar values. We did not account for depletion of gas-phase HN0 3 
or HC1 in the calculation of the equilibrium compositions; it should not be significant 
for background aerosol sulfate loadings in the temperature range considered in this 
work. 16 For ease of discussion, we will hereafter refer to each solution by its equilibrium 
temperature and by whether it contains nitric acid [ i.e ., 197 K (no HN0 3 )]. 

In Fig. 4 we present loss and production curves for C10N0 2 and Cl 2 , respectively, 
for the 198 K (no HN0 3 ) solution; the two curves yield similar first-order rate con- 
stants. The reaction probability values (y) for the two nitric acid scenarios are listed in 
Tables 2 and 3, and are the result of an average of at least three separate measurements. 
Through the use of bulk solutions prepared with specific compositions, we have effec- 
tively decoupled the atmospheric equilibrium temperature from the laboratory tem- 
perature; the former has a much larger effect on the reaction probabilities. Because of 
the tendency of the supercooled solutions to freeze, it was necessary to measure the 
reaction probabilities at laboratory temperatures which are higher than the atmospheric 
equilibrium temperatures. As shown in Table 4, the variation in the reaction probability 
for the 197 K (no HN0 3 ) solution over the laboratory temperature range 203-233 K 
was quite small, of the order of the experimental uncertainty. The 195 K (HNO a ) solu- 
tion lies well within the nitric acid trihydrate (NAT) stability region of the phase 
diagram, 8 and we were unable to cool the solution below 223 K without freezing. All 
other reaction probability measurements were carried out at 213 K and usually also at 


Table 2 Reaction probabilities for the ternary H 2 S0 4 -H 2 0-HC1 system measured at 
203 K, as a function of atmospheric equilibrium temperature 


T(eq)/K weight fraction H 2 S0 4 weight fraction HC1 reaction probability 


201 

0.552 

3.00 x 10“ 6 

0.010 

199 

0.520 

1.60 x 10“ 5 

0.022 

198 

0.510 

2.40 x 10“ 5 

0.038 

197 

0.490 

6.00 x 10“ 5 

0.061 

196 

0.470 

1.40 x 10“ 4 

0.083 

195 

0.450 

3.44 x 10“ 4 

0.12 

194 

0.430 

8.24 x 10“ 4 

0.30 


Table 3 Reaction probabilities for the quaternary H 2 S0 4 -HN0 3 -H 2 0-HC1 system 
measured at 203 K, as a function of atmospheric equilibrium temperature 

T(eq)/K 

weight fraction 

h 2 so 4 

weight fraction 

hno 3 

weight fraction 
HC1 

reaction 

probability 

199 


0.022 

1.38 x 10 5 

0.021 

198 


0.035 

2.70 x 10“ 5 

0.026 

197 


0.061 

6.18 x 10“ 5 

0.058 

196 

0.362 

0.126 

1.76 x 10“ 4 

0.070 

195 


0.283 

4.88 x 10“ 4 

0.110 
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Fig. 4 #, Loss (CIONO,) (Ac = 315 s *) and ■, production (Cl 2 ) {k = 322 s *) curves for the 
HC1 + C10N0 2 reaction on the 198 K (no HN0 3 ) solution measured at a laboratory temperature 

of 213 K 


203 K, so that the extrapolation to the atmospheric temperature was not extreme. 
Hanson and Ravishankara 5 had already noted that the reaction probability is nearly 
independent of temperature for solutions of a given composition; they attributed this 
behaviour to a fortuitous compensation of effects. 

In Fig. 5, the results from our reaction probability measurements are presented and 
compared with values estimated from the results reported by Hanson and 
Ravishankara 5 for ternary (H 2 S0 4 -H 2 0-HC1) compositions which closely correspond 
to our 196, 198 and 201 K solutions. It may be seen that our reaction probabilities for 
solutions without HN0 3 are slightly lower than those of Hanson and Ravishankara. 5 
Fig. 5 also shows that the incorporation of HN0 3 into the sulfate solutions lowers the 
reaction probability relative to solutions without HN0 3 only to a small extent, despite 
the presence of HNO a levels as high as 28.3%; the difference is of the order of the 
experimental uncertainty, although it appears to be statistically significant. Zhang et al. 6 
also measured reaction probabilities without HN0 3 at a temperature of 198 K, for one 
solution containing ca. 4% HN0 3 , and found no change in the C10N0 2 + HCl reac- 
tion probability. 

In Fig. 6 the reaction probabilities from Tables 2 and 3 are instead plotted as a 
function of HCl content. The line in the figure is a least-squares fit; its slope indicates 
that y is directly proportional to the square root of the HCl concentration. Hanson and 
Ravishankara 5 measured y over a wide HCl partial pressure (p HC1 ) range: as expected, at 
very low p HCI the C10N0 2 decay rate is dominated by reaction (1), and hence is nearly 
independent of p HC , . For intermediate and high p HC1 values y was found to be pro- 
portional to (p H ci) 1/2 an d to p HC1 , respectively; Hanson and Ravishankara 5 interpreted 
these results in terms of a bulk and a surface component to y. The measurements of 
Zhang et al 6 were carried out at relatively high HCl concentrations, and their data 
indicate close to linear proportionality between y and p HCI . Considering the validity of 
Henry’s law for these conditions, the same functional forms should apply to [HCl] : our 
measurements were carried out at intermediate [HCl] values, yielding a square-root 
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0.14 


0.12 


0.10 


Lo 

2 0.08 
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O 

a. 

| 0.06 
o 

CO 

£ 

0.04 


0.02 


0 

194 195 196 197 198 199 200 201 202 

7-(eq)/K 

Fig. 5 Reaction probability measurements for HC1 + C10N0 2 as a function of atmospheric tem- 
perature T(eq), assuming 5 ppm by volume H 2 0 and 2 ppb by volume HC1 at 100 mbar. □, Our 
data, for 10 ppb by volume HN0 3 ; ■, our data, assuming no HN0 3 ; ▲, data estimated from the 
results of Hanson and Ravishankara 5 (no HN0 3 ). 



proportionality, in agreement with the results obtained by Hanson and Ravishankara 5 
for similar [HC1] values. 

In an attempt to ascertain a dependence of y on the H 2 S0 4 content, we chose a fixed 
HC1 concentration of 2.3 x 10" 3 mol dm" 3 [corresponding to the moderately efficient 
reaction probability of the 197 K (no HN0 3 ) solution] and varied the H 2 S0 4 content 
from 40 to 60%. This data is listed in Table 4; it is clear that the reaction probability 



Fig. 6 Reaction probability measurements for HC1 + C10N0 2 as a function of HC1 weight frac- 
tion. □ and ■, experimental data for solutions with and without HN0 3 , respectively. The line is 

a fit through the data; its slope is 0.545. 
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Table 4 Reaction probabilities for the H 2 S0 4 -H 2 0-HC1 
system as a function of temperature and sulfuric acid content, 
measured at a fixed HC1 concentration of 2.3 x 10“ 3 mol dm - 3 


T(exp)/K 

weight fraction H 2 S0 4 

reaction probability 

203 

0.49 

0.060 

213 

0.49 

0.061 

223 

0.49 

0.050 

233 

0.49 

0.049 

213 

0.40 

0.057 

213 

0.45 

0.046 

213 

0.51 

0.061 

213 

0.55 

0.059 

213 

0.60 

0.054 


does not change significantly over the stratospherically relevant H 2 S0 4 content range 
given in the table. This result, in combination with the laboratory temperature- 
dependence studies, suggests that the atmospheric temperature (and thus the composi- 
tion of the aerosols) affects the C10N0 2 + HC1 reaction probability predominantly 
through its effect on the HC1 solubility. Although the uptake of HN0 3 serves to increase 
the solubility of HC1, it may be seen (Tables 2 and 3) that reaction probabilities are 
already very high at the atmospheric temperatures where this effect becomes important. 

Conclusions 

The Henry’s Law solubility constants for HC1 reported in this work for the ternary 
H 2 S0 4 -H 2 0-HC1 and quaternary H 2 S0 4 -HN0 3 -H 2 0-HC1 systems have been shown 
to agree with the thermodynamic model of Luo et a/. 13 for selected atmospherically 
relevant compositions and temperatures. Furthermore, our new results for the ternary 
system are in agreement with the previous measurements of Hanson and Ravishankara 4 
and Williams and Golden, 11 while our earlier results 12 deviate significantly from these 
consensus findings. 

Our results show that the uptake of substantial levels of nitric acid (as much as 28%) 
does not appreciably affect the reaction probabilities for the C10N0 2 + HC1 process on 
sulfuric acid surfaces. The results also show that the reaction probabilities are not sensi- 
tive to temperature for solutions with a given composition. The controlling parameter is 
the concentration of HC1 in the solution, which for a given HC1 partial pressure depends 
in turn on temperature; in the lower stratosphere this concentration increases rapidly 
below ca. 200 K. 

Since it is now clear that the C10N0 2 4- HC1 reaction proceeds efficiently on super- 
cooled liquid sulfuric acid surfaces with compositions relevant to the cold stratosphere, 
including the possibility of extensive HN0 3 uptake, it is not necessary to invoke the 
presence of solid PSC particles to explain the high chlorine activation rates associated 
with large scale ozone depletion. As suggested earlier, the key factor is the stratospheric 
temperature; 8 as it approaches the HN0 3 -H 2 0 dew point, the composition of the 
liquid aerosols changes rapidly, 22 with the consequent increase in HC1 solubility. If the 
aerosols are solid, it is the approach to the ice frost point that matters; chlorine activa- 
tion occurs readily on the ice-like surfaces. We note, however, that a change in phase for 
the stratospheric aerosols might have other important consequences: the surface area 
available for reaction can change with freezing and, more importantly, the mechanism 
for efficient denitrification might involve nucleation of a small fraction of the liquid 
particles; the solid ones would then selectively grow and fall out of the stratosphere. 
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